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1 .  INTRODUCTION 


If  the  conductivity  and  dielectric  constant  of  two  distinct  media 
are  known,  the  amount  of  electromagnetic  energy  which  will  be  reflected 
at  their  interface  can  be  related  to  the  incident  energy  through  the 
Fresnel  equations.  On  a  large  scale,  the  earth-air  boundary  is  an 
example  of  such  an  interface.  In  comparison  to  the  earth,  the  electri¬ 
cal  properties  of  the  air  can  be  considered  relatively  independent  of 
regional  variations.  Therefore,  obtaining  typical  electrical  properties 
of  the  soil  in  a  particular  region  becomes  a  primary  objective  if  the 
reflection  coefficients  are  desired. 

With  regional  values,  the  reflection  of  electromagnetic  energy,  such 
as  that  encountered  after  a  high-altitude  nuclear  burst,  can  be  calcu¬ 
lated.  Without  regional  values,  however,  the  spectrum  of  energy 
reflected  may  be  inadequately  predicted.  Any  method  of  approximating 
conductivity  and  dielectric  constant  will  produce  an  error  if  it  ex¬ 
cludes  regional  variables  such  as  structure  or  mineral  content,  or  if  it 
excludes  frequency  and  moisture  content.  To  determine  whether  this 
error  should  cause  significant  concern,  the  actual  values  for  conductiv¬ 
ity  and  dielectric  constant  can  be  measured  and  compared  with  predicted 
values.  In  this  way,  the  reflected  energy  spectrum  can  be  more  accu¬ 
rately  predicted. 


2.  SOIL  PARAMETERS 

In  order  to  determine  the  energy  of  a  pulse  wave  reflected  from  the 
surface  of  the  earth,  the  electrical  properties  of  this  soil/air  bound¬ 
ary  must  be  known.  The  Fresnel  equations,  which  describe  the  reflection 
of  electromagnetic  energy,  require  that  the  conductivity,  permittivity, 
and  permeability  of  both  the  air  and  the  soil  be  known.  Under  normal 
conditions,  and  in  comparison  to  soil,  the  electrical  properties  of  air 
are  approximately  equal  to  those  of  free  space.  Therefore,  determining 
the  electrical  properties  of  the  earth  becomes  a  primary  concern. 

Several  studies  have  been  conducted  concerning  the  electrical  prop¬ 
erties  of  the  soil.  The  results  of  these  studies  state  that  the  conduc¬ 
tivity  (o)  and  permittivity  (e)  can  be  predicted  at  a  specified 
frequency,  if  the  moisture  content  is  known.  The  permeability  (p), 
which  is  a  function  of  the  ferrous  mineral  content,  deviates  slightly 
from  that  of  free  space  in  most  regions. 

Scott1  measured  numerous  samples  over  the  frequency  range  of  102  to 
106  Hz.  After  averaging  his  data,  he  produced  a  set  of  curves  for 


d.  Scott,  Electrical  and  Magnetic  Properties  of  Rock  and  Soil, 
U.S.  Geological  Survey,  Note  18  (1966), 
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conductivity  and  dielectric  constant  as  functions  of  frequency.  These 
curves,  which  have  become  known  as  Scott's  universal  curves,  have  been 
curve  fit  with  a  second-order  surface  fit  and  are  mathematically  repre¬ 
sented  in  equations  ( 1 )  and  ( 2 ) . 


Scott's  curve  fit  for  conductivity  is 

K  =  -0.604  +  1.640W  -  0.062F  +  0.062W2  -  0.070FW  +  0.021F2;  (1) 

the  curve  fit  for  relative  dielectric  constant  is 

D  =  4.905  +  1.308W  -  0.971F  +  0.111W2  -  0.160FW  +  0.059F2,  (2) 

where 

K  =  log1Q  of  conductivity  (mmho/m), 

D  =  log1Q  of  dielectric  constant  (e/eQ), 

W  =  log1Q  of  water  content  (percent  by  volume),  and 
F  =  log10  of  frequency  (Hz). 

Since  Scott's  work  encompasses  a  large  variety  of  laboratory- 
measured  soil  samples  (samples  which  correlated  well  with  field  measure¬ 
ments),  equations  (1)  and  (2)  give  an  average  of  conductivity  and  di¬ 
electric  constant  over  various  regions. 

C.  L.  longmire2  developed  a  time-domain  representation  called  tht 
universal  RC  network  model  (fig.  1).  The  basic  assumption  of  this  model 
is  that  the  soil  can  be  regarded  electrically  as  an  equivalent  network 
of  resistors  and  cap?citors.  Longmire  made  no  laboratory  measurements, 
but  by  using  Scott's  curves  and  data  that  he  received  by  private 
consultation  with  Wilkenfeld,2  he  extended  the  expected  validity  of  his 
model  over  the  frequency  range  of  10°  to  3  x  1012  Hz. 


Figure  1.  Universal  network  model. 


2C.  L.  Longmire,  A  Universal  Impedance  for  Soils,  Defense  Nuclear 
Agency,  DNA  3788T  (1975), 


Longmire's  model  states  that 
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where 

C  is  capacitance, 
a  is  a  constant, 

R  is  resistance, 

Rq  is  the  resistance  of  the  circuit  that  is  asymptotically  ap¬ 
proached  as  frequency  is  decreased, 
o  is  conductivity, 
f  is  frequency,  and 

the  subscript  «°  indicates  a  value  which  is  asymptotically  approached 
as  frequency  is  increased. 


The  relative  dielectric  constant  of  the  soil  is 
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and  the  conductivity  is 
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(3) 


N  (f/fn)2 

0=0+  2 ire  l  a  f  - - -  (mho/m)  ,  (4) 

0  °  n=1  n  n  1  +  (f/fn)2 

where  fn  accounts  for  moisture  content.  Table  1  gives  coefficients  for 
universal  soil  as  derived  by  Longmire's  model. 

TABLE  1.  COEFFICIENT  an  FOR  UNIVERSAL  SOIL 
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To  obtain  a  standard  for  evaluating  the  error  which  might  be  in¬ 
curred  by  applying  the  results  of  this  network  model  in  a  region  which 
does  not  follow  the  typical  soil  parameters,  laboratory  measurements 
were  made  of  er  and  o  for  five  samples  of  soil  taken  from  the  Harry 
Diamond  Laboratories,  Woodbridge  Research  Facility,  Woodbridge,  VA.  The 
National  Bureau  of  Standards  conducted  these  measurements  over  the 
frequency  range  of  5  x  1 06  to  6  x  1 08  Hz.  Following  Scott's  work,  a 
second-order  surface  fit  was  made  for  these  data  as  given  in  equations 
( 5 )  and  ( 6 ) : 

D  =  5.886  -  1.045F  -  2.055W  +  0.077FW  +  0.056F2  +  1 . 1 80W2  ,  (5) 

K  =  -3.233  +  0.6453F  -  0.448W  -  0.105FW  -  0.0113F2  +  0.8876W2  .  (6) 


Equations  (5)  and  (6)  are  henceforth  referred  to  as  the  varying  regional 
parameter  equations. 

The  averaged  values  at  5-,  10-,  15-,  20-,  and  25-percent  moisture 
content  for  conductivity  and  dielectric  constant  are  shown  in  figures  2 
and  3.  Also  shown  are  the  values  predicted  by  the  Longmire  universal 
network  model  for  10-percent  moisture  content  (eq  (3)  and  (4)). 

As  can  be  seen  from  figure  3,  averaged  values  for  conductivity  are 
considerably  below  the  values  predicted  by  the  network  model.  Since 
these  data  were  obtained  in  a  region  of  higher  than  average  rainfall, 
and  taken  from  the  uppermost  soil,  it  is  believed  that  their  values 
indicate  only  that  the  surface  soil  possesses  a  low  conductivity. 
Several  other  field  measurements  were  taken  near  the  laboratory, 
measurements  which  indicate  that  the  conductivity  of  the  soil  increases 
sharply  a  few  feet  below  the  surface  of  the  earth.  Unfortunately,  the 
samples  were  too  few  to  determine  where  the  sharp  increase  took  place 
nor  what  conductivity  the  soil  might  possess  in  the  lower  strata. 

When  either  conductivity  or  dielectric  constant  changes  from  one 
homogeneous  medium  to  another,  as  is  the  case  at  the  Woodbridge  Research 
Facility,  a  second  medium  for  reflection  is  provided.  Because  these 
boundaries  are  close  to  each  other,  the  energy  which  is  transmitted 
through  the  primary  boundary  will  undergo  a  second  reflection  at  the 
secondary  boundary.  The  wave  which  will  ultimately  be  seen  by  an  ob¬ 
server  will  be  a  combination  of  an  incident  wave,  a  primary  reflected 
wave,  and  a  series  of  multiply  reflected  waves.  This  problem  is 
discussed  briefly  after  the  reflection  from  the  primary  boundary  is 
discussed. 
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3.  REFLECTION  COEFFICIENTS  FROM  SOILS 


3.1  Single  Reflection 


Soils,  like  most  common  materials,  fall  in  the  category  of 
lossy  dielectrics.  Fresnel’s  equations  (as  discussed  in  many  electro¬ 
magnetics  texts3)  describe  the  reflection/transmission  process  of  elec¬ 
tric  and  magnetic  components  of  an  electromagnetic  field  across  the 
boundaries  of  two  distinct  media.  The  following  equations  relate  to 
figures  4  and  5. 

For  the  E  field  normal  to  the  plane  of  incidence. 


U2k1 (COS  02  +  C°S  V 
E!  "  cos  02  +  M1k2  COS  91  E°  * 


(7) 


3J.  A.  Stratton,  Elect romgnetic  Theory,  McGraw-Hill  Book  Co,  (1941). 


For  the  H  field  normal  to  the  plane  of  incidence 


(cos  02  +  cos  0Q) 

H!  ~  U^2  cos  @2  +  y2^1  cos  01  H0  ' 


y1^2  cos  0Q  -  y2^1  cos  ©1 
H2  "  y1k2  cos  0Q  +  y2k1  cos  0,|  Ho  ' 


where  0  is  angle  of  incidence,  and  k^  is  the  complex 
constant  defined  by 


k?  =  » 

where  (i)  =  2irf. 

Through  Snell’s  law 

sin  0q  =  sin  02  ,  and 

k2  sin  0q  =  k1  sin  01  . 

The  reflection  coefficients  can  be  written 

y2k1  cos  0Q  -  -  **  sin2  0Q 

y2k1  cos  0Q  +  y, yjk \  -  k^  sin2  0Q 
for  E  normal  to  the  plane  of  incidence,  and 

*Sk2  COS  90  ~  W2k1>/k2  ~  k?  Sin  90 

W1k2  COS  90  +  y2k1>/k2  _  k?  Sin  90 
for  H  normal  to  the  plane  of  incidence. 


(9) 


(10) 


propagation 


By  considering  a  typical  air/earth  interface,  these  formulas 
can  be  further  simplified.  Medium  2  (characterized  by  ©2 ,  P2'  e2^  *-s 
that  of  the  earth,  and  medium  1  (characterized  by  o1 ,  ,  e1 )  is  that  of 

air;  since 


“l  “  “2  “  ^0 


£1  =  £0  ' 


0-,  -  o  , 


the  reflection  coefficients  can  be  simplified  to 


cos  8 


cos  0Q  + 


v)  ~si 

0  -S‘ 


2  O 

sin  6, 


•  2  a 
sin  0, 


V)C°S  9°  ~  j  V")  ~  31 

/  °>  \  FI  ~M 

U  -  j  - )  cos  0„  +  /(e  -  j  - )  -  si 

\ r  £q V  0  V\r  0  / 


sin2  0, 


•  2  A 

sin  0, 


where  er  =  £i/£q  and  subscripts  H  and  V  refer  to  horizontal  and  vertical 
components. 

Figures  6,  7,  8,  and  9  show  the  magnitudes  and  phases  of 
reflection  coefficients  RH  and  Ry  using  varying  regional  parameter 
equations.  The  angle  of  incidence  is  60°. 

The  incident  and  reflected  waves  combine  to  form  what  is 
generally  referred  to  as  the  total  wave.  These  equations  have  been 
developed  in  other  texts*4  and  are  restated  here: 

It 

M.  A.  Messier ,  The  Effects  of  Ground  Reflection  on  Observed  EMP 
Waveforms,  Defense  Nuclear  Agency,  DNA  3370T  (1974). 


E. 


E 


E_ 


Eo  L1  *R»e  J 
Eo  E  +  V'3  <h/>>] 
Eo[’  +  V'J(3/X)] 


COS  (j)  , 


sin  $  cos  9  , 


sin  $  cos  9  , 


(13) 

(14) 

(15) 


where 

X  =  wavelength, 

h  =  distance  between  point  of  reflection  and  observer, 

4>  =  polarization  of  E  vector,  parallel  to  vector  n, 

=  radial  wave,  and 
Eq  =  incident  wave. 

These  equations  apply  to  the  geometry  of  figure  10  and  assume  that  the 
distance  between  the  source  of  the  pulse  and  either  the  observer  or  the 
earth  is  much  greater  than  the  distance  between  the  earth  and  the  ob¬ 
server. 


Figure  6.  Magnitude  of  reflection  coefficients  (RH)  using  varying- 
parameter  model  at  5-,  10-,  15-,  and  20-percent  moisture. 
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Figure  7.  Pha'>  of  reflection  coefficients  (RH)  using  varying-parameter 


model  at  5-,  10-,  l5-#  and  20-percent  moisture. 
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Figure  8.  Magnitude  of  reflection  coefficients  (Ry)  using  varying- 
parameter  model  at  5-,  10-,  15-,  and  20-percent  moisture. 
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model  at  5-,  10-,  15-,  and  20-percent  moisture. 
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Figure  10.  Singly  reflected  wave. 
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3.2  Multiple  Reflection 


If  there  is  a  marked  increase  in  soil  conductivity  at  some 
ith,  is  is  the  case  for  most  of  the  eastern  seaboard,  5  a  transmitted 
i.  ill  undergo  a  primary  reflection  and  a  series  of  multiple 

reflec  ions.  Figure  11  depicts  this  process.  The  total  wave  seen  by  an 
observer  just  above  the  earth  will  be  composed  of  an  incident  wave,  a 
primary  reflected  wave,  and  a  series  of  multiply  reflected  waves.6 

The  multiply  reflected  pulse  is  a  geometric  series: 


E  =.  e  [  1  +  f o  +  TT^O'*  e^2 

r  oL  v.v  h 

+  TT'p'2pMej2A2  +  rr'p'V2  ej3A2  ...)]  (16> 

(for  zero  polarization) 
where 

Er  =  resultant  E  field, 
p  =  reflection  coefficient, 

T  =  transmission  coefficient, 

A1  and  A2  account  for  phase  delay  and  attenuation  in  medium  1  and 
medium  2,  respectively,  and 

primes  indicate  the  number  of  times  the  wave  has  been  reflected. 

PRIMARY 

REFLECTED  MULTIPLY  REFLECTED 
WAVE  WAVES 


12  3  4 


Figure  11.  Multiply  reflected  wave. 


5E.  D.  Sunde,  Earth  Conduction  Effects  in  Transmission  Systems,  D.  Van 
Nostrand  Co,  Inc .  (1968). 

6tf.  V.  Klein,  Optics,  T.  Wiley  and  Sons  Inc .  (1970). 


Equation  (16)  can  be  reduced  to  the  closed-form  equivalent: 


FREQUENCY  (Hz) 

Figure  13.  Magnitude  of  multiple  reflection  coefficient  (rh)  for  a  two- 
layer  earth.  Varying-parameter  model  is  used  for  upper  layer  (1  m 
thick),  and  network  model  is  used  for  lower  layer  (infinite  half  plane). 


FREQUENCY  (Hz) 

Figure  14.  Magnitude  of  multiple  reflection  coefficient  (RH)  for  a  two 
layer  earth.  Varying-parameter  model  is  used  for  upper  layer  (2  m 
thick),  and  network  model  is  used  for  lower  layer  (infinite  half  plane) 


FREQUENCY  (Hz) 

Figure  15.  Magnitude  of  multiple  reflection  coefficient  (Ry)  for  a  two- 
layer  earth.  Varying-parameter  model  is  used  for  upper  layer  (0.5  m 
thick),  and  network  model  is  used  for  lower  layer  (infinite  half  plane). 


FREQUENCY  (Hz) 

Figure  16.  Magnitude  of  multiple  reflection  coefficient  (Ry)  for  a  two 
layer  earth.  Varying-parameter  model  is  used  for  upper  layer  (1  m 
thick),  and  network  model  is  used  for  lower  layer  (infinite  half  plane) 
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Figure  17.  Magnitude  of  multiple  reflection  coefficient  (Ry)  for  a  two- 
layer  earth.  Varying-parameter  model  is  used  for  upper  layer  (2  m 
thick),  and  network  model  is  used  for  lower  layer  (infinite  half  plane). 


4.  ERROR  ANALYSIS 

A  convenient  way  to  represent  the  total  energy  on  a  per  ohm  basis  is 
by  integrating  the  magnitude  of  the  spectral  density  over  the  range  of 
frequencies  of  interests 


r-..  i 


E  “  17  /L  |EC«) I2  dw  . 


The  energy  difference  between  two  real  spectra  is 


7  Iq  |  E(u>)  |  ^  -  |  E(c*>)  |  ^  dw  . 


(26) 


(27) 


Here,  E(w)  represents  the  spectrum  of  the  total  E  field  of  the  signal 
used  as  the  basis  for  analysis,  and  E(«)  is  the  approximating  signal. 
The  percentage  of  error  (Pe)  is 
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(28) 


/"  { |  E<  to)  |  2  -  |i(0J)|2}  dw 
p  =  - - - x  100% 

®  /"  |  E (to )  |  2  do) 

Since  it  is  not  possible  to  numerically  integrate  over  the  entire  spec¬ 
trum,  the  analysis  will  include  the  extrapolated  spectrum  of  105  to  109 
Hz  in  which  the  variable  parameter  model  is  valid.  Also,  to  eliminate 
the  confusion  or  cancellation  due  to  phase  shift  caused  by  the  delay  of 
the  reflected  wave,  the  observer  is  placed  an  incremental  distance  above 
the  earth's  surface.  Here  also,  the  analysis  is  done  for  zero  polar¬ 
ization  (i.e.,  Rv  is  not  used).  Table  2  lists  the  percentages  of  error 
between  the  total  energies  of  waves  above  soil  modeled  by  the  variable- 
parameter  model  and  soil  modeled  by  the  network  model  with  the  variable- 
parameter  model  as  the  standard. 

The  values  show  that  the  spectra  E(w)  and  E (to)  are  displaced  from 
each  other  and  will  probably  give  values  for  Pg  which  will  be  no  less 
than  50  percent. 

Another  model  which  is  often  used  for  calculating  propagation  con¬ 
stants  and  reflection  coefficients  is  the  constant-parameter  model. 
This  model,  as  the  name  implies,  assumes  that  the  conductivity  and 
dielectric  constant  are  constant  over  the  entire  spectrum  of  frequencies 
and  range  of  possible  moisture  contents.  Typical  values  for  these 
parameters  are  0  =  7  mmho/m  and  er  *  15  for  the  Woodbridge  Research 
Facility.  The  values  for  the  percentage  of  error  are  listed  in  table  3. 

These  values  indicate  that  by  judicious  selection  of  0  and  er  the 
percentage  of  error  can  be  reduced.  However,  some  error  is  inevitable 
because  the  spectrum  predicted  by  the  constant-parameter  model  will  have 
a  sharper  roll-off,  as  demonstrated  in  the  comparison  of  the  three 
models  given  in  figure  18. 


TABLE  2.  PERCENTAGE  OF  ERROR  FOR  VARIABLE-PARAMETER  MODEL 
(STANDARD)  AND  NETWORK  MODEL  FOR  VARIOUS  MOISTURE  CONTENTS 


Moisture  content  (%) 

Percentage  of  error  (%) 

5 

59.0 

10 

-75.2 

15 

-80.4 

20 

-82.7 
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TABLE  3.  PERCENTAGE  OF  ERROR  FOR  VARIABLE- 
PARAMETER  MODEL  (STANDARD)  AND  CONSTANT- 
PARAMETER  MODEL  FOR  VARIOUS  MOISTURE  CONTENTS 


Moisture  content  (%) 

Percentage  of  error  (%) 

5 

68.9 

10 

63.0 

15 

54.4 

20 

43.8 

FREQUENCY  (Hz) 

Figure  18.  Comparison  between  reflection  coefficient  as  predicted  by 
varying  regional  parameter  model,  network  model  (both  at  10-percent 
moisture),  and  constant-parameter  model  (conductivity  =  7  mmho/m, 
dielectric  constant  =  15). 


5.  CONCLUSION 


The  percentage  of  error  for  the  energy  difference  between  two 
signals  was  minimized  by  using  the  constant-parameter  model.  This, 
however,  does  not  mean  that  the  constant-parameter  model  is  the  best 
possible  method.  Ihe  spectra  of  the  reflected  wave  had  a  sharper  roll¬ 
off  than  either  the  variable-parameter  model  or  the  network  model.  This 
means  that  there  will  be  some  minimum  error  which  cannot  be  eliminated. 

If  there  were  data  between  102  and  106  Hz,  the  network  model  could 
be  altered  to  more  closely  approximate  the  true  nature  of  the  soil.  A 
convenient  method  for  obtaining  the  soil  parameter  at  any  frequency  and 
moisture  content  would  then  exist.  It  could  conceivably  minimize  the 
percentage  error  below  10  percent. 

However,  the  problem  generated  by  the  multiple  reflection  from  a 
stratified  medium  would  still  remain.  The  only  solution  for  this  would 
be  to  follow  the  stated  equations  for  multiple  reflections  after  the 
lower  stratification  has  been  characterized  by  extensive  field  measure¬ 
ments. 
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